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Helical Luttinger liquids, appearing at the edge of two-dimensional topological insulators, repre¬ 
sent a new paradigm of one-dimensional systems, where peculiar quantum phenomena can be in¬ 
vestigated. Motivated by recent experiments on charge fractionalization, we propose a setup based 
on helical Luttinger liquids that allows to time-resolve, in addition to charge fractionalization, also 
spin-charge separation and pure spin fractionalization. This is due to the combined presence of 
spin-momentum locking and interactions. We show that electric time-resolved measurements can 
reveal both charge and spin properties, avoiding the need of magnetic materials. Although challeng¬ 
ing, the proposed setup could be achieved with nowadays technologies, promoting helical liquids as 
interesting playgrounds to explore the effects of interactions in one dimension. 


I. INTRODUCTION 

Many body physics in low dimensions is a very ac¬ 
tive field in condensed matter 1 . Here one-dimensional 
(ID) interacting systems play an important roltP®. since 
different fascinating phenomena have been predicted 1 -®® 
and observed^®! The breakdown of the Fermi liquid 
paradigm in ID led to the introduction of the Tomonaga- 
Luttinger liquid theorjtf^UJU. This theory predicts pecu¬ 
liar phenomena, such as the spin-charge separation, i.e. 
interacting electrons split up into charge and spin col¬ 
lective excitations with different velocities 1 2 . This has 
been measured in semiconductor quantum wired^, carbon 
nanotube^ and quasi-ID organic conductor^. Another 
intriguing aspect is charge fractionalization: an electron 
injected into an interacting region breaks up in collec¬ 
tive excitations which carry fractional charge^®®®! 
Although many theoretical works have concentrated on 
this, experimental evidences have been elusive for a long 
time. A first observation was achieved by St einberg 
et al. via momentum-resolved spectroscop y 10 ® ! More 
recently, Kamata et al. have reported the first direct 
measurement of fractional excitations by means of time- 
resolved transport measurements in a Hall bar setup 1 !®! 
In this work, we propose a setup in which, in addition 
to charge fractionalization and spin-charge separation, 
also pure spin fractionalization can be observed. This 
is not expected to occur in ordinary ID liquids, due 
to spin-rotation invarianc^^. Our proposal is based on 
the edge states of two-dimensional (2D) topological in- 
sulatonPH®, whose properties are described in terms 
of counter-propagating channels characterized by spin- 
momentum locking. In the presence of interactions, these 
are described as helical Luttinger liquids (hLLjPUH, 
whose first experimental manifestations have been re¬ 
cently observed®! Here, charge and spin sectors are 
connected and the presence of interactions implies the 
simultaneous fractionalization of both degrees of free¬ 
dom 30 36 . However, evidence of pure spin fractionaliza¬ 
tion, i.e. fractional spinons with no charge, have not 


yet been observed. Indeed, the presence of measuring 
contacts prevents the observation of fractionalization via 
standard dc transport measurement:! 19 ®! Therefore, we 
propose time-resolved detection sc heme , whose degree of 
control has been greatly improved 11 3 ®aH. However, in 
hLL spin-momentum locking^ 4 ®! prevents the indepen¬ 
dent observation of their fractionalization 1 . To circum¬ 
vent this limitation, we consider a gate electrode to un¬ 
balance charge and spin velocities, allowing to separately 
identify fractionalized holons and spinons. Our proposal 
relies only on a time-resolved electrical scheme, where 
also spin properties can be probed by means of conven¬ 
tional charge current measurements. The proposed setup 
would allow to detect three fingerprints of interactions in 
the ID world (spin-charge separation and their fraction¬ 
alization) in a single experiment. These results provide 
an important step forward in the understanding of the 
nature of interactions in one dimension. 


II. MODEL 

We consider the setup in Fig. [IJa)-(b), with helical 
edge modes of a 2D topological insulator described by 

Hoix) = tj;\{x){-ihv F d x )tp^(x) + tl>l(x)(ihv F d x )'ipi(x). 

(!) 

Short-range interactions are present, given by 
'H {L \ x ) = i ) ^2 P°( x ) [54110*0 At (a) + g 2 ±{x)p- a (x)] 

a=t,I 

( 2 ) 

with pc, =: 'i/jJ.'i p a : the electron particle density of each 
channel. By negatively polarizing the central electrode 
the edge states are forced to change their profile. Here we 
assume that the gate electrode is wide enough to prevent 
tunneling, but narrow enough to allow short-range inter¬ 
actions between the edge states on the two arms. This is 
indeed the case in the experiment performed by Kamata 
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FIG. 1. (Color online) (a) Spin up (red) and spin down (blue) 
electrons of a single helical edge state counter propagate. 
Electron wave-packets can be pumped into the system from 
the injector and probed by the detector. An additional gate 
electrode (yellow) changes the edge state profile, (b) Model: 
the white regions correspond to the Fermi liquid leads, the 
gray ones to the edge states away from the gate electrode and 
the yellow ones to the edge states around the gate electrode, 
(c) Folding procedure applied to (b), that maps the infinite 
2-channel (helical) liquid into the semi-infinite 4-channel 
(spinful) one. 


et al. in Ref. EH Therefore the presence of the central 
electrode induces additional non-local interactions 


= )- ^2 Pa(-x) [g 2 \\ {x)p a (x) + g 4 ±(x)p- a (x)\ . 
<r=t,4 

( 3 ) 

The spatial dependence of the interaction parameters 
is 9a\\(2±)(x) = g 4 \\ (2_L) Q(T/2 - M) and g 4 _L(2||) (x) = 
5 , 4 _ l ( 2 ||) 0((/2 — |a;|). The edges are connected to Fermi- 
liquid metallic contacts, modeled as ID non-interacting 
system^®. Note that a rigorous microscopic calculation 
of the dependence of the strength of interactions on the 
geometric gate parameters would involve the full numer¬ 
ical solution of a bidimensional Poisson problem posed 
by the whole setup, which is beyond the scope of the 
present article. Nevertheless, modelling the effect of the 
gate voltage on the edge state physics through short- 
range Coulomb interactions in the form of Eq. [3] has 
been shown to be a very good effective description, cor¬ 
roborated by experimental results 11 . 

Interactions are treated within bosonization 1 ^. Introduc¬ 
ing the scalar fields 4> a related to the electron density as 
p^[/\.{x) = /ii x ) the Hamiltonian is diagonal¬ 

ized in terms of the fields 


/ <t>c{x) 
9 c (x) 
<t>s{x) 
V 0 s {x) 


-11 1 -1 \ / <f> t(:r) \ 

111 1 ^t(-z) (4) 

-11-1 1 Mx) ’ 1 j 

i i —i —i / V / 



sLL-region 

0 < x < 1/2 

hLL-region 

1/2 < x < L/2 

FL-region 
L/2 < x 

K^(x) 

K c 

K h 

1 

K {s) (x) 

K s 

1/AT 

1 

v^(x) 

Vc 

Vh 

Vf 

v( s \x) 

Vs 

Vh 

Vf 


TABLE I. Luttinger parameters and velocities appearing in 
the inhomogeneous Luttinger liquid model Eqs. , cor¬ 

responding to the setup of Fig. [T] 


satisfying [d x $\(x), 9\>{x')\ = iS x ,x'S(x - x'). The 
bosonized version of H = j ^ dx [Hq(x) +i~L^ L \x) + 
W l \x)\ can be written as H = Ylx=c s with 


H x = 


dx 


M\x) 


d x (j)\{x) 


KW{ x) 


K (x \x ) d x O x (x) 


1 2 


(5) 

The transformations Eq. ([4]) fictitiously double the de¬ 
grees of freedom of the system, whose correct number is 
restored by halving the domain of integration. The in¬ 
finite 2-channel liquid with non-local interactions, Fig. 
0b), is mapped onto the semi-infinite 4-channel liquid 
with local interactional, Fig. [ljc). Equations (4| de¬ 
fine charge and spin fields of the spinful-like 4-channel 
liquid, related to its charge (units e) and spin (units 
fr/2) densities and currents as p x = x/2jird x (f)\ and 
jx = — x/2/irdt^x- Equation (ij) represents the Hamil¬ 
tonian for a spinful Luttinger liquid (sLL). Charge and 
spin sector separates, characterized by velocities 


^ (A) (aO = xfsJ( 1 + g 4 ,x{x)f - glx(x) (6) 

and Luttinger parameters 


kM(: r) = 


'1 + g 4 ,x(x) - 92,xi x ) 

1 + 94,\{ x ) + 92 ,x(x)’ 


(7) 


with M*) = a "' ( 1L 8 T l,) and *,<*) = 

We refer to 1/2 < x < L/2 as the hLL-region, and sLL- 
region the one for 0 < x < 1/2, see Fig. 0c). As 
shown in Table [TJ inhomogeneous interactions generate 
space-dependent parameters. This inhomogeneous Lut¬ 
tinger liquid model has been applied t o exp lain conduc¬ 
tance quantization in interacting wires l 5 l 6 l 15 l and c harge 
fractionalization in quantum Hall device^- 11 20 The 

setup of Fig. 0 a ) represents a particular inhomogeneous 
Luttinger liquid, where Fermi liquid (FL), hLL and sLL 
behaviors coexist. At the FL-hLL interface (hLL|FL) 
at x = L/2 and at the hLL-sLL interface (sLL|hLL) at 
x = 1/2 , the inhomogeneous interactions affect the prop¬ 
agation of excitations. 
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III. TIME EVOLUTION 

We study space-time evolution of electron wave- 
packets injected from the left contact (injector) and de¬ 
tected by the right one (detector). Injection can be 
achieved by e.g. applying a voltage pulse to the injec¬ 
tion gate, while quantum point con tacts ca n be exploited 
to perform time-resolved detectior l 11 l 39 l 43 l ljH Crucially, 
the concept of electron wave-packet makes sense in the 
FL-region only: once entered the interacting region, the 
electron wave-packet is decomposed into charge and spin 
collective excitations, governed by the equations of mo- 
tiorpEEH 


dUx(x,t) = v^(x)K^(x)d x 

obtained by recalling 
v < ' X \x)K < ' X \x)d x 6\(x,t) and 

(x) 


V ( X ) 

KW(x) 


d x (j>x{x,t ) 


( 8 ) 


d t d\(x,t) 


KW(x) t) ■ From Eq. ( 8 ) one can deter¬ 

mine the dynamical evolution of <p\(x 1 t) = ((f>\(x,t)}, by 
imposing the continuity of <f>\(x,t) and d x <f>\(x, t) 

in x = 1/2, L/2, together with d t (f>\(0,t) = 0, the 
latter imposing total reflectioiP^ at x = 0 (see Fig. 
[ljc)). The equation of motion can be mapped into a 
scattering problem for an incident plasmon mode, as 
described in Appendix [A] Space and time derivatives 
of 4>\(x,t) are related to the average densities p\(x,t) 
and currents j\(x,t) respectively. To determine their 


dynamical evolution, the initial (t = 0 ) density p^ (x) 
and current j^\x) profiles, determined by the injection 
process, must be specified. Because of helicity, only 
spin-up electrons can propagate from the injector. 


Therefore we consider p~c\x) = pi°\x) = p®( x) 


and jc°\x) = js V> (x) = VfP^°\x), with 

p(°)(x) = Vexp[(x — Xi) 2 /2cr 2 ]/(v / 27rcr) assumed to 
have a Gaussian distribution, X{ > L/2 being the 
injection point. N is the number of injected elec¬ 
trons, with total charge Q and spin S = hN/2. The 
parameter a can describe the finite duration of the 
injection process; it is related to the full-width-at-half- 
maximum St of the time distribution of the injected 
wave-packets as a ~ 0A2v F St. By solving Eq. the 
dynamical evolution of charge and spin sectors can be 
determined (see Appendix Aj) . This can be written as 
p\{x,t) = p^ nc \x,t) + p ( ^{x,t), with p { ™ c) (x,t) the 
= wave-packet incoming from the injector and p^ e ^ l \x,t) 


AO), 


the wave-packet reflected to the FL-region (see Fig. 
[ljc) ). The latter encodes information about scattering 
phenomena and fractionalizations, and is related to the 
detected current (units e) as 


jdet(x d , t) = ~ p^ e//) (x d , t) + p { s refl) (x d , t ) 


2 L‘ 


(9) 


with x d > L/2 the detection point. Crucially, due to he¬ 
licity, information about both charge and, more impor¬ 
tantly, spin can be extracted from the detected current. 
By recalling Eqs. ( A28|A29 ) one has 


p^ efl \x d ,t) = rp {0 \-Xd + v F t + L) +7 E C' njm p (0) (-Xd + VF t + L - (n + 1) — (L - /) - m—l) 

nti 0 V V h V C ) 


( 10 ) 


— 


1 -K h 
1 + K h 


K c - I< h 

I<c + K h 


n+m+1 


min{m,n-(-l} 

(1 - S m , 0 ) E 

P=1 


m—l 

p- 1 


n + 1 
P 


~4K h K c 


(k c - K h y 


"m, 0 


( 11 ) 


r 


with r = 


7 = 


4 K h 

' ( 1 +K h ) 2 


The spin density evolu¬ 


tion pi re f l \x d , t) can be obtained from Eqs. 00 with 
the substitutions 


(v c ,K c ,K h ) (v„K„l/K h ). 


( 12 ) 

Note that in the non-interacting case py e ^ l \x d ,t) = 
pi re ^ l \x d ,t) = p^(—x d + v F t): the injected wave-packet 
simply propagate through the system with constant ve¬ 
locity vf, and does not undergo neither fractionalizations 
nor spin-charge separation (see inset of Fig. [2]). 

The first term in the r.h.s. of Eq. (10) corresponds to the 


injected wave-packet being reflected at hLL|FL. Because 
of helicity, the reflection of a fraction rQ of the injected 


charge Q is accompanied by a reflected fraction —rS of 
the total incoming spin S. Spin fractionalization is due to 
the helical nature of the edge states, which, in the pres¬ 
ence of interactions, effectively breaks the spin-rotation 
symmetry. The component transmitted across hLL|FL 
generates the second term of Eq. (10), where an infinite 
number of terms appears. Once entered the interacting 
region, the transmitted wave-packet can undergo multi¬ 
ple reflections between the interfaces before going back to 
the FL-region. Terms proportional to C n , m correspond 
to n reflections at hLL|FL toward hLL-region and m re¬ 
flections at x = 0 (see Appendix [A|. 

Note that for long times / 0 °° dtji re ^ l \x d ,t) = Q 
and / 0 °° dtji refl \x d ,t) = S, with j^ efl \x d , t) = 
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v F Px efl \x d ,t). Then, from Eq. |d|, / 0 °° dtj det (x d ,t) = 
Q: in the dc limit, the injected wave-packet is trans¬ 
mitted to the detector, in agreement with the universal 
quantization of dc conductance^. 


IV. RESULTS AND DISCUSSION 

The detected current as a function of time jdet(xd,t) 
is shown, with realistic parameters, in Fig. [2j together 
with charge and spin currents of the 4-channel liquid. 
Different signals are revealed at different times, reflecting 



FIG. 2. (Color online) Detected, charge and spin (multiplied 
by 2e/h) currents as a function of time. The injected wave- 
packet consists of N = 10 spin-up-polarized electrons, whose 
time distribution has a width-at-half-maximum 5t as 40 ps. 
Parameters are Xi = Xd = 12 p m, L = 18 /im, l = 10 p m. 
The Luttinger parameter^ 49 ! are v c = vf/K c , vh = vf/Kk, 
v s = vfK s , vf = 3 x 10 4 m/s and K h = 0.6, K c = 0.7, 
K s = 1.1. The inset shows the detected current without 
interactions, where neither fractionalizations nor spin-charge 
separation are present. 


signatures of fractionalization. We focus on the signals 
a-d, that are captured by the short-time contributions of 
Eq. (10). For the charge sector they read 


p£ ef i\x,t) « rp^\—x + VFt + L ) 

+ 7Co,oP (0) + v F t + L - VF —y 


+ 7 Co,i/^ 0) ( —x + v F t + L — v F - 


L — l 
Vh 


V F - 


(13) 


and similarly for the spin sector, by using Eq. (12). The 


first charge and spin components, first member in the 
r.h.s. of Eq. (13), corresponds to the injected wave- 


packet being reflected before entering the interacting re¬ 
gion, and thus are not probed by the detector (signals a). 
The components transmitted across hLL|FL propagate 
in the hLL-region as charge and spin density waves. At 
sLL|hLL they are partially reflected back, second member 
in the r.h.s. of Eq. jlil), with charge +jc, ) Q 


and spin ( i ( k^k s + i ) i leading to a finite detected 
current (peak b). However, this signal does not allow to 
independently identify charge and spin fractionalization: 
the current Eq. (|9]) is a superposition of both fractional 
charge and spin. 

The transmitted components across sLL|hLL, third mem¬ 
ber in the r.h.s. of Eq. (13), produce signals c and 
d in Fig. [2j The corresponding physical processes are 
shown in Fig. [3] Once the wave-packet enters the sLL- 



FIG. 3. (Color online) Physical processes corresponding to 
signals c and d in Fig. [2] Charge and spin excitations are rep¬ 
resented in green and orange respectively, (top) The charge 
and spin components of the injected wave-packet are frac¬ 
tionalized when entering the interacting region. Here they 
propagate with the same velocity vh reaching the gate elec¬ 
trode, where they are transmitted and further fractionalized. 
Around the electrode charge and spin separate ( v c > v s ). 
(bottom) The charge (spin) component returns back and is 
split into two wave-packets, represented by dashed green (or¬ 
ange) rectangles, with opposite spin (charge), one propagat¬ 
ing toward the detector and the other one toward the injector. 
The components propagating toward the detector are finally 
probed, with a time shift At, as electric signals in the time- 
resolved current (solid rectangles). 


region, its charge and spin components separate, with 
charge density waves propagating faster than spin ones 
(v c > v s ), as shown in Fig. [3](top). Therefore the charge 
excitation returns back to sLL|hLL before the spin one. 
Here the charge (spin) wave-packet is split into two wave- 
packets with opposite spin (charge), marked by dashed 
green (orange) rectangles in Fig. [^bottom). The compo¬ 
nents propagating toward the detector are finally probed 
as electric signals c and d in Fig. [2j with a time delay 
At = 1 due to the different time of flights of pure 

charge and pure spin excitations. This time shift rep¬ 
resents a manifestation of spin-charge separation. The 
peak marked with c in Fig. [2] is due to fractionaliza¬ 
tion of pure charge excitations (holons), while the peak 
marked with d is related to fractionalization of pure spin 
(spinons). The other signals observed in Fig. [2] are 
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due to several reflections inside the interacting region. 
Note that, for different interaction parameters, the pre¬ 
cise structure of the current can change, but the signals 
reflecting pure charge and pure spin fractionalization can 
be identified in any case. 

Therefore, we come to the main result of our work. 
The proposed setup allows to observe the landmarks 
of interactions in one dimension by means of a unique 
electric time-resolved detection scheme: fractionalization 
of holons and spinons can be identified, the time delay 
in their detection providing manifestation of spin-charge 
separation. This is particularly fascinating for the possi¬ 
bility to electrically probe pure spin fractionalization. 


V. EXPERIMENTAL FEASIBILITY 

Key point of the proposal is the ability to time-resolve 
the charge and spin fractionalization, signals c and d in 
Fig. 0 This implies At > St , that is, the spin-charge 
time shift At greater than the width of the time distribu¬ 
tion of the wave-packets St. The ability to inject electron 
wave-packets with St of the order of tens of piRA together 
with realistic values of the Fermi velocity (vf ~ 3 x 10 4 
m/s for the edge states in InAs/GaSb quantum wells), 
gives a lower bound on the system size of approxima- 
tively 1 /iin. An upper bound is given by the inelastic 
mean free path Z; n . Ballistic transport has been observed 
up to few /im-long samples, but further comprehension 
of the scattering mechanisms 50 51 would allow to increase 
the coherence length up to tens of /rniSUSsH 57 ! The plot 
of Fig. [2] corresponds to ~ 10 /im-long setups, and 
the charge and spin signals, although slightly overlap¬ 
ping, can be resolved. Therefore, the proposed measure¬ 
ments can be likely performed with current technologies, 
and helical edge states reasonably represent an promis¬ 
ing playground to explore the effects of interactions in 
one dimension. 
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Appendix A: Time evolution of charge and spin 
density 

Here we concentrate on the time evolution of the av¬ 
erage quantities <px(x,t) = ((f>x(x,t)) which satisfy the 


equations of motion 


dUx (x,t) = vW(x)K^(x)d x 


'/ X ^(x) 


d x 4>\{x,t ) 


KW{x) 

(Al) 

The parameters v^ x \x) and K^ x \x) are piecewise 
constants, being discontinuous at the junctions be¬ 
tween different regions. To evaluate the densities 
px(x,t) = x/2/'Kd x (f)\[x 1 1) and currents jx{x,t) = 

— x/2/ndt4>\(x, t ) it is necessary to evaluate the space- 
time evolution of <j>\ (x,t). This can be done solving Eq. 

(Al) by imposing that <j>\{x,t) and ^.^?\ d x (j>\{x,t) are 
continuous functions and th at t he current = 0. 

Due to the linearity in Eq. </>x{x,t) can be decom¬ 

posed into bosonic modes 


<l>\,q{x, t) = e lVFqt Sx, q {x) 
with momentum q and energy hvFq as 

p-\-oo 


/ -too 

dq a\ t q (f>X,q{x, t), 

-oo 


(A2) 


(A3) 


the coefficients ax, q being related to the distribution of 
the injected wave-packet, as will be shown in the follow¬ 
ing. Therefore, the problem of solving the equation of 
motion for <f>x(x,t) is mapped into a scat terin g problem 
for the incident mode <j>x,q(x,t) in Eq. (A2), in agree¬ 
ment with the generalized plasmon scattering approach. 
In the following we focus on the charge sector only, and 
omit the index A = c for notational convenience (for 
example s q (x) = Sa, 9 (x)); the spin sector with A = s 
can be obtained with the substitutions (■ v c , K c , ) —> 
(v s ,K s ,l/K h ), as reported in Table I of the main text. 
For the charge sector 


j(x) = 


K(x) = 


V F 

x > L/2 


Vh 

1/2 < x < L/2 

(A4) 

v c 

0<x<l/2 


1 

x > L/2 


K h 

1/2 < x < L/2 . 

(A5) 

I<c 

Q<x<l/2 



Therefore to solve the scattering problem one requires 
that 


r e~ iqx + R q e iqx x > L/2 

s q (x) = l A q e~' iq ^ x + B q e iq ^ x l/2<x<L/2 , 
{c q e~ iq ^ x + D q e iq ^ x 0 < x < 1/2 

(A6) 

and k(x) d x s q (x) are continuous functions at x = 
1/2, L/2, together with s g (0) = 0 which guarantees 
j(0,t) = 0. With these constrains it is possible to deter¬ 
mine R q , A q , B q , C q and D q . In particular, we are inter¬ 
ested to evaluate R q , which is related to the components 
of the injected wave-packet entering the detector (in the 
four-channel model of Fig. 1(c) in the main text, only 
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reflected components of the injected wave-packet can be 
prob ed b y the detector). Once R q is known, we can use 
Eq. (A3) to get, for x > L/ 2, 


p{x,t) = \ - 


J~ r°° dq a q e~ iqVFt d x (e~' qx + R q e iqx ) 
y ^ J —oo 


= p^{x,t) + p^ refl \x,t), 


(A7) 


where the incoming 


/ o /*+°° 

- J dq a q qe- iq ^ x+vpt) (A8) 

and reflected 

P {refl) (: x , t) = i^J + dq a q qR q e iq{x ~ Vpt) (A9) 


components have been identified. At t = 0 the incom¬ 
ing wave-packet is specified by the injection process, i.e. 
p(* nc )(x, 0) = p(°\x), which is non-vanishing close to the 
injection point xi in the FL region. By comparing Eq. 
(A8) with the Fourier transform of the injected wave- 
packet 


/ +oo 

-OO 


— /j(°) e~ iqx 
2n Pq ’ 


(A10) 


one finds a q = y j^ K - 2q P'q' 1 an d from Eq. (A9) 

p {refl \x,t) = - J + °° J p^R q e iq ^- VFt) . (All) 


By means of Eq. (All) it is possible to determine the 
evolution of the reflected component from the knowledge 
of p q >> , related to the shape of the injected wave-packet, 
and from R q , related to the solution of the scattering 
problem, that is explicitly determined in the following. 


1. Reflection coefficient R„ 


Solving the scattering problem is, in principle, straight¬ 
forward: the five coefficients appearing Eq. (A6) can 


indeed be determined from the five boundaries condi¬ 
tions. However this approach is not convenient from a 
physical point of view. Indeed, it does not allow to dis¬ 
criminate the different physical processes occurring in the 
system, due to multiple reflections of the injected wave- 
packet. To overcome this limitation and be able to follow 
the space-time evolution of the collective excitations, we 
adopt a different strategy. We c onsider all possible paths 
that a bosonic mode Eq. ( |A2| ), once entered the inter¬ 
acting region x < A/2, can follow before being reflected 
back in the region x > L/2. A possible example of one 
of these paths is shown in Fig. E] We call Ji the hLL|FL 
junction at x = A/2, J\ the sLLjhLL junction at x = 1/2 
and Jq the end of the whole (folded) system at x = 0. 


m — 3 

Jo 

" <z 


€ 


C < 


j i 


n = 3 

J 2 

-©-<- 


J 


J 




L_ 
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X 


FIG. 4. (Color online) An example of a possible path con¬ 
tributing to the reflection coefficient R q . After the injected 
wave packet is transmitted across J 2 to the interacting region 
(upper purple circle) and before it is reflected back to the 
FL-region x > L/2 (lower purple circle), it undergoes mul¬ 
tiple reflections. In particular, the specific path illustrated 
here is characterized by n = 3 reflections at J 2 (blue dots) 
and m = 3 reflections at Jo (red dots). There are p = 2 
disconnected subpaths inside the sLL-region x < 1/2. There¬ 
fore there are 2p = 4 transmissions across Ji (green circles), 
m — p = 1 reflection from Ji toward Jo and n + 1 — p = 2 
reflections from Ji toward J 2 , in agreement with Eq. (|A13|). 


The first contribution to the reflection coefficient R q is 
obviously due to the first reflection from J 2 : the injected 
wave-packet does not enter the interacting region. We de¬ 
note this contribution as r q the arrow (—►) denoting 
incoming wave-packet incident from the right (left). As 
for the other coefficients that will be introduced in the 
following, we refer to the next section for its explicit eval¬ 
uation. The other contributions are due to the fraction of 
the injected wave-packet transmitted across J 2 into the 
interacting region x < L/2. Here the fractionalized wave 
packet can undergo several reflections from Jo, Ji and J 2 
before going back to x > L/2 , where it can be probed 
by the detector. All these contributions pick up a factor 
tq^t q ^/l, i.e. the product of the transmission coefficient 
across J 2 when entering (t q J /l) and exiting (t q J ^}) the in¬ 
teracting region (see upper and lower purple circles in 
Fig. [^respectively). Note that, once entered the inter¬ 
acting region x < L/2 , the wave packet arrives at J\ and, 
just before exiting the interacting region, it comes from 
Ji again (see Fig. |4]) . In order to consider all the con¬ 
tributions to R q we conveniently focus on the paths that 
start and arrive at J\. 

Let n be the number of reflections at J 2 toward J\ and 
m the number of reflections at Jo (toward J\ obviously). 
Both n and m can assume natural values. For the specific 
case shown in Fig. [4] n = 3 and m = 3. The reflection 
coefficient from J 2 toward J\ is r q J ?}, while with r q J /l 
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we denote the reflection coefficient from Jq. The contri¬ 
butions to R q due to several reflections is thus given by 

and (r q ^j . Therefore the global reflection 
coefficient can be written as 

+ OO 

£ (+?>}(+i>) a,,. 

n,m=0 

, (A12) 

The coefficient P n ,m takes into account the contribution 
of all the different paths with fixed m and n. Indeed, 
assigning (n, m) does not completely specify the path. 
Moreover, also transmission and reflections at J\ con¬ 
tribute to the coefficient of a specific path, and are en¬ 
coded into P n ,m- Let us see this point in more detail. 

If m = 0 there are no reflections from Jo and P Ut o is given 

by ( 'n h ll 1 , i.e. the n + 1 reflections at J\ toward J 2 . 

The calculation is more complicated if to > 0. It is use¬ 


^ min{m,n+l} 

Pn,m = S m , 0 ( r ^) + (1 “ £m,o) 

P= 1 



To determine the global reflection coefficient Eq. (A12), 
we need to find the explicit expressions of the reflection 
and transmission coefficients through each barrier. This 
is done in the next Section. 


2. Reflection and transmission coefficients at the 
barriers 


Here we explicitly determine the coefficients r q J /l , r^a 
and tq J a with a =<—. —> and i = 1,2. Let us consider a 
junction located at x = xq, which separates regions A 
(for x < xq, characterized by Ka and va) and B (for 
x > Xq characterized by Kb and vb)- This situations 
applies both to Ji (xq = 1/2 ) and J 2 (xo = L/ 2). The 
scattering states are 


u 

u 


iq^-x , —iq^-x 

e V A +r q ^e v a 

X < Xo 

t q ^e iq ^ x ’ 

X > Xo 

t q M~ iq ^ X 

X < Xo 

— iq-^-x . iq—^~ x 

e + r 9i< _e 

X > Xo 


corresponding to modes incident from A ( f q 
B respectively. The continuity of f q ^(x) 


u(x 


K(a 


rf q ^(x) implies 


and 

and 


iq V F(A[-^)xo 

Ka + I<b 

I<b - Ka i2 q ^x 0 
-e 11 a u . 

K a + K b 


(A14a) 

(Al4b) 


ful to introduce the integer number p , that counts the 
number of disconnected subpaths inside the sLL-region 
(.x < 1/2). For example, p = 2 in Fig. [3] The number p is 
bounded as 1 < p < min{m, n + 1}. If there are p discon¬ 
nected subpaths in the sLL-region, then there are p trans¬ 
mission across Ji from the hLL-region to the sLL-region 

j , p transmissions across J\ from the sLL-region to 
the hLL-region , ( m—p ) reflections at J\ toward 


J 0 and (n + TO +1) — 2p— (m — p) = n+1 — p 

f ,j 

reflections at Ji toward J 2 ( r q .ll 1 


Finally, the grouping of the m and n reflections into p 
subpaths is not unique, so that additional combinatory 
factors must be kept into account. The final expression 
for P n ,m reads 


("; I )( r ™)" +1 - p (a)”- r (++<«)'. <ai3) 


r 


In the opposite case, for scattering from the right the 
coefficients are 


tq ^ - K A + I< n 


2K A ^iqv F (/j-^)x o 


Ka — K b _* 2 q ^£- 


K a + K b 


i2q-^-xo 


(Al5a) 

(Al5b) 


Using the expressions in Eqs. (A14i-b) and (A15r-b) it 
is easy to show that 


Ml) _ 1 -iqL 

9 K h + 1 

r (+) = 

9 1 + K h 

dJ 2 ) _ 2K h 1)/2 

K h + 1 

(W = 2 e l9L{ A7- l)/2 

9 ’-*■ K h + 1 

Mi) = Kc ~ Kh 
9 K c + I< h 

r ( + ) = Kh ~ K c r i q l |g 


K c + K h 
2 K r 


Mi) — 

9 K c + K h 


O y V v c vu" 


= 

9 K c + K) 


2Kh iq u VE_3LE)/ 2 

-e K v c v h J/ . 


(A16) 

(A17) 

(A18) 

(A19) 

(A20) 

(A21) 

(A22) 

(A23) 


In addition, focusing on the scattering state incoming on 
Jo and requiring that the current vanishes at x = 0 one 




















finds 


r (Jo) — 
q,<- 


- 1 . 


(A24) 


Substituting Eqs. ( A16|A24 ) into Eqs. (A12) and (A13)) 
we obtain 


~~ K 


K h - 1 


-iqL 


4 K h 


JqL^-1) 


(K h + l ) 1 2 3 4 5 6 7 8 


oo 

E 

n,m =0 


(-D r 


1 - K h 

1 + Kb 


iqL-^-n 71 

e v h P n 


(A25) 


and 

Pn,m — ^m, 0 


K c - K h \ ” +1 _ igZ iE( n+1 ) 
K c + K h ) 


, C \ -ifci|f(ra+l) ikl^im) (K c — K h \ + + J /to - 1\ fn + l\ / -AK c R h \ 

+ (1 ~ <W)e ' e •• UtkJ '- 11 g U-ijl „ KnwipJ 


min{ra,n+l} 


Using straightforward algebra the global reflection coefficient can be expressed in compact form as 


Rn = ^—^e- iqL + exp 


9 K h + 1' 


(K h + l ) 2 


iq ( (L — l) — (n + 1) + l — m — L 
Vh v c 


(A26) 


(A27) 


with 


Cn 7 m 


1 -K h \ n (K c -K h 


1 + K h 


Kr + Kb 


n+m+1 


min{ra,n+l} 

(1 - o) 

P= 1 


?n — 1\ /n + 1 
P-1 


-K< h K c 


i v 


(K c - 


+ S', 


m,0 


Finally, by restoring the index A = c for the charge sector, using Eq. (All) we find 


OO / \ 

pi refl) (x d ,t) = rp (0 \-x d + v F t + L) + 7 V] Cn m p (0) ( + A - (n + 1)^:|L - () - m—l ) 

V v h v c ) 


(A28) 


(A29) 


with r = gpfg, 7 = ~ ( 1 4.^)2 ■ The spin density evolution p s (x d ,t) can be obtained from Eqs. ( A28JA29 ) with the 
substitutions (v c ,K c ,K h ) -t (v s , K s , 1/K h ). 
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